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ABSTRACT 

We investigate the acceleration of charged particles (both electrons and protons) at collisionless 
shocks predicted to exist in the vicinity of solar flares. The existe nce of standing termination shock s 
has been examined by flare models and numerical simulations (e.g.. lShibata et al. 1995; Forbcs fl986ft . 
We study electron energization by numerically integrating the equations of motion of a large number 
of test-particle electrons in the time-dependent two-dimensional electric and magnetic fields generated 
from hybrid simulations (kinetic ions and fluid electron) using parameters typical of the solar flare 
plasma environment. The shock is produced by injecting plasma flow toward a rigid piston. Large- 
scale magnetic fluctuations - known to exist in plasmas and known to have important effects on the 
nonthermal electron acceleration at shocks - are also included in our simulations. For the parameters 
characteristic of the flaring region, our calculations suggest that the termination shock formed in the 
reconnection outflow region (above post-flare loops) could accelerate electrons to a kinetic energy 
of a few MeV within 100 ion cyclotron periods, which is of the order of a millisecond. Given a 
sufficient turbulence amplitude level (SB 2 /B 2 ~ 0.3), about 10% of thermal test-particle electrons are 
accelerated to more than 15 keV. We find that protons are also accelerated, but not to as high energy 
in the available time and the energy spectra are considerably steeper than that of the electrons for the 
parameters used in our simulations. Our results are qualitatively consistent with the observed hard 
X-ray emissions in solar flares. 

Subject headings: Sun: flares - acceleration of particles - shock waves - turbulence 



1. INTRODUCTION 

Solar flares are obser yed to be strong so urces of ener- 
getic charged particles (jAschwandenll200"2ft . The release 
of magnetic energy by mag netic reconnection is thought 
to be the driving process (|Masuda et al.l[l99l . While 
several mechanisms have been proposed to explain the 
acceleration of ch arged particles in flares (see review by 
iMiller et al.lll997l and references therein) , there is still no 
general consensus and this remains an unsolved problem. 
Recent hard X-ray observations of the nonthermal elec- 
tron bremsstrahlung emission by Reuven Ra maty High 
Ener gy Solar Spectroscopic Imager (RHESSI; iLin et all 
12002ft have provided more details of electron accelera- 
tion in solar flares. The observations indicate a large 
fraction of released energy resides in high-energy elec- 
trons accelerated soon after the flare. Hard X-ray sources 
above the top of magnetic loops have bee n detected (e.g., 
IMasuda et"aUll994tlKrucker et all 12010ft . providing im- 
portant clues to the acceleration process. For example, 
the lo op-top source recently reported by IKrucker et al.l 
(|2010ft shows a large number of electrons (> 5 x 10 35 ) 
are accelerated to more than 16 keV and the highest en- 
ergy reaches ~ MeV. Accelerated ions have also been 
observed in solar flares and tend to correlate with en- 
ergetic electrons (jShih et al.l 120091 ) . Since the observed 
hard X-ray source requires a very efficient acceleration, 
explaining how such a large number of electrons (prob- 
ably also ions) are accelerated to high energy poses a 
challenge to theoretical astrophysics. 

The existence of fast shocks in the reconnection out- 
flow region has been predicted in flare models (e.g., 

guofan@lpl.arizona.edu 



Shibata et al.l 119951) and numerical s i mulations (e.g., 
Forbes! 119861 . [1988: Sh iota et al.l l2Tjol IWorkman et al.l 



20111 ). Using MHD numerical simulations. iForbesI (|1986ft 



studied the formation of a standing termination shock 
when a high-speed jet driven by reconnection encoun- 
ters a closed magnetic loop. The geometry of the flare 
termination shock can be represented by Figure [TJ The 
high-speed jet created in the reconnection outflow region 
collides with the top of the magnetic loop and produces 
a fast-mode, standing termination shock. The resulting 
flare termination shock has a unit normal to its surface 
that points nearly perpendicular to the magnetic field. 
This is a perpendicular shock (i.e., the angle between 
the upstream magnet i c field and shock normal vector 
&Bn = 90°). IForbesI (|1986D predicts the existence of 
this shock with a compression ratio of 2.0 and an up- 
stream Mach number as high as 2.3. A recent study by 
IWorkman et al.l (|2011l ) shows similar results. The obser- 
vational evid ence of the exi s tence of flare shocks has been 
presented bv lAurass et al.l (|2002l ). 

Fast-mode, collisionless shocks are known to 
be efficient accelerator s of charged particles 
(jBlandford fc Eichlerl 11987ft. The theory of di ffusive 
shock acceleration (e.g.. lAxford et al.l 119771; iBelll 11978ft 
predicts shock-accelerated particles have a power-law 
distribution, which is often seen in solar energetic particle 
events. Energetic electrons are frequently observed to be 
associated with collisionless shocks with 0r„ > 45° (e.g.. 
Fan et al.l 119641: lAnderson et al.l 11979) : iTsurutani fc Linl 
19851; iSimnett et all 120051: lOka et all 120061) . Electron 



acceleration by obliq ue shocks has been considered by 
a number of author s (iHolman fc Pesseslll983l; IWul 11984 ; 



iLeroy fe Mangeneyl 119841 ; iKrauss-Varban et all 



1989; 



2 



Guo and Giacalone 



Krauss-Varban fc Wul Jl98l IGuo fe Giacalonel [20101: 
Riguelme fc Spitkovskvl 120111: iMatsukivo et al.l 1201 ID . 
In the scatter-free limit, the acceleration of electrons at 
a planar quasi-perpendicular shock is often referred to 
as fast-Fermi acceleration or shock drift acceleration. 
It can be demonstrated that these two mechanisms are 
the same process in two di fferent frames of reference 
(|Krauss-Varban fc Wul |1989() . Shock drift acceleration 
at a flare termination shock in the reconnection outflow 
region has been considered by a number of authors (e.g., 
iMann et afl 120091: IWarmuth et aD l2009f l . However it is 
commonly known that in this process both the fraction 
and attainabl e energy of the acce lerated particles are 
limited (e.g.. iBall fc Melrose! [200lh . The acceleration 
and scattering of low-rigidity particles, such as electrons 
is thought to be a problem. This is commonly referred to 
as "injection problem" . For low-rigidity ions, the injec- 
tion problem is better understood. It is usually thought 
that the ion injection at quasi-parallel shocks is easier 
tha n that at quasi - perpen dicular shock (see discussions 
by iLembege et al.l (120041) and more re cent numerical 
simula tion studies bv lYang "et~aT1 ([2009] ) and IWu et al.l 
<|2009f ) ) . However, numerical simulations that consider 
large-scale pre-existing mag netic turbulence suggests 
there is no injection problem (|Giacalo nc 2005ajb|). There 
has been some works on facilit ating particle accel e ration 
at flare termination shocks. iTsuneta fc Naitd (|1998ft 
considered electron heating by a slow-shock p air as 
a pre-energization process. iSomov fc Kosugil (|1997l ) 
considered the role of plasma heating and collapsing 
magnetic trap at reconnecting magnetic field lines. 

One of the recent advances in our understanding of 
electron acceleration at shocks is the consideration of 
nonuniform effects such as shock ripples a n d mag netic 
turbulence. The s imulations by iBurgess! (|2006l f and 
lUmeda et al.l (|2009l ) showed that small-scale shock rip- 
ples can be important in sc attering the electrons and 
facilitating the acceleration. iSavoini fc Lembegel (|2010T ) 
studied the motions of electrons in the vicinity of a shock 
using two-dimensional full particle simulations. They 
found the effects of shock reformation and shock rippling 
are important for non-adiab atic electron heating at col - 
lisionless shocks. Recently, Uokipii fc Giacalone] (|2007f ) 
proposed a mechanism to solve the injection problem 
which does not require strong pitch-angle scattering from 
small-scale fluctuations. The fast-moving electrons can 
move along meandering magnetic field lines and travel 
back and forth across the shock front, and, therefore 
gain energy from the difference between the upstream 
and downstream flow velocities. Using self-consistent hy- 
brid simulati on combined with test-p article simulation 
for electrons, IGuo fc Giacalone] (|2010T l have found effi- 
cient electron acceleration at perpendicular shocks mov- 
ing through a plasma containing large-scale pre-existing 
upstream magnetic turbulence. The turbulent magnetic 
field leads to field-line meandering which allows the elec- 
trons to get accelerated at the shock front multiple times, 
like in the Jokipii & Giacalone picture. The shock sur- 
face becomes irregular on a variety of spatial scales from 
small-scale ripples due to ion-scale plasma instabilities 
(|Lowe fc Burge ss 2003), to large-scale structures caused 
by the interaction bet ween the shock and upstream 
large-scale tur bulence (jGiacalone fc Neugebauerl 120081 : 
ILu et al.l 12009ft . The rippled surface of the shock front 



also contributes to the acceleration by mirroring elec- 
trons between the ripples. The observational evidence of 
these lar ge-scale ripples have been shown by a number of 
authors (INeugebauer fc Giacalonel2005UBale et al.lll999t 
iPulupa fc Bald 120081 : iKoval fc Szabol l2010f ). These re- 
sults, along with t he previous work on acceleration of ions 
(|Giacalone 2005a. b), suggest that large-scale turbulence 
has an important effect on the acceleration of both elec- 
tron and ions at shocks, which is consistent with the cor- 
relation between i ons and elec t rons in solar energetic par - 
ticle events (e.g.. iCliverl 120091 : iHaggertv fc Roelojl2009f) . 
This mechanism especially favors perpendicular shocks, 
particularly for electrons, and could explain the observed 
high-energy particles in solar flares as well. 

In this paper we present results from a combination of 
hybrid simulations and test-particle simulations to study 
the electron and proton energization at a shock with pa- 
rameters similar to the flare termination shock inferred 
by previous works in the existence of upstream magnetic 
fluctuations. We consider the nonlinear modification of a 
planar shock front by upstream Alfvenic fluctuations and 
its effect on electron acceleration. Although the plasma 
waves and turbulence in the reconnection outflow region 
could be considerably different from our simplified model, 
some intrinsic characteristics of this interaction, such as 
the braiding of magnetic field lines and shock rippling 
are still preserved. Our study suggests the acceleration 
of electrons at a flare termination shock is rapid and ef- 
ficient. A large fraction of the initial thermal electrons 
is accelerated to hundreds of keV and even reaches MeV 
energies in a very short time. This indicates collisionless 
shocks may play an important role in particle accelera- 
tion in solar flares. In section 2 we describe the numerical 
method used in this paper. Section 3 discusses the simu- 
lation results. In section 4 we discuss the implication of 
our results and draw conclusions. 

2. NUMERICAL METHOD 

We use a combination of two-dimensional hybrid simu- 
lations and test-particle electron simulations to study the 
energization of electrons at the flare terminat ion shock. 
In the hybrid simulation model (Winskc 1985), the ions 
are treated kinetically and the background thermal elec- 
trons are usually treated as a massless, charge neutraliz- 
ing fluid. The electrons are solved by using electron fluid 
momentum equation (generalized Ohm's law). This ap- 
proach is well suited to resolve the kinetic physics of ions 
and has been successfully applied to study the physics 
of collisionless shocks. We consider a two-dimensional 
Cartesian geometry in which all physical quantities are 
functions of two spatial variables x and z. All vectors 
have components in three directions (x, y, z). The size of 
the simulation domain L x x L z for each case is listed in 
Table [1] The flare termination shock is modeled by in- 
jecting plasma continuously from one end (x = 0) of the 
simulation box and colliding with a rigid wall at the other 
end (x = L x ). A simplified one-dimensional fluctuating 
magnetic field B(z, t) = <5B(z, t) + B\qz is present in the 
simulation box at t = and also injected continuously 
at the x = boundary with the plasma flow during the 
simulation, where Bio is the averaged upstream magnetic 
field. The fluctuating component (5B(z,i) consists of an 
equal-partition of right- and left-hand circular polarized, 
forward and backward parallel-propagating plane Alfvcn 
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waves. The amplitude of the fluctuations at the wave 
number k is determined from a Kolmogorov power spec- 
trum: 



where L c is the coherence length. The total variance, 
which is the integral of P over all k, in each case is listed 
in Table 1. In all the cases, we take the turbulence coher- 
ence length to be equal to the box size in the z direction 
L c = L z . We expect large-scale magnetic turbulence to 
exist in the reconnection outflow plasma. These fluctua- 
tions can be triggered by reconnection, foot -point motion 
or other processes fsee lKrucker et all . 120081 . for a detailed 
discussion). Although the simplified form of magnetic- 
field fluctuations we use may not be representative of the 
reconnection outflow plasma, any turbulence with large 
variances should allow strong field-line wandering, which 
is essential in our particle acceleration model. 

In order to produce a low Mach number shock, as pre- 
dicted by other numerical simulations, the inflow Alfven 
Mach number is taken to be Mao = 1-0. After reflec- 
tion at the right boundary, this produces a shock with 
an average Mach number of about 2.0 in the shock frame, 
consistent with the flare te rmination shock predicted by 
previ ous MHD simulations dForbc s 1986; Wor kman et al.l 
[201 If). The grid sizes are Ax = Az = 0.5c/w P i and the 
time step is taken to be At = O.Olfi^ , where c/u> P i is the 
ion inertial length and fi" 1 is the ion cyclotron period. 
The ion and electron plasma betas, /3j and /3 e , respec- 
tively, are both taken to be 0.03 and the ratio between 
light speed and upstream Alfven speed c/Vao = 410, 
which roughly corresponds to an initial situation with 
proton and electron temperatures T,; = T e = 2 x 10 6 A', 
number density n = 8 x 10 9 cm~ 3 and magnetic field 
Bp = 30G, similar to constraints from observations 
(|Krucker et al.l 12010ft . For these parameters, the aver- 
age upstream speed in shock frame of reference (also the 
jet outflow speed) is about 1460km/ s. This speed is mea- 
surable because the outflow plasma from reconnection is 
moving at the Alfven speed. The estimate from obser- 
vations (|Tsuneta et al.l 11997ft is roughly consistent with 
this value. 

Since in the hybrid model the electrons are treated as 
a massless fluid, the test-particle simulations are needed 
to study their kinetic behavior. In the second step, we 
integrate the relativistic equations of motion for an en- 
semble of test-particle electrons in the two-dimensional 
time-dependent electric and magnetic fields obtained in 
the hybrid simulations. These test-particle electrons are 
treated as a different part from the electron fluid in the 
hybrid simulations. We use a second-order spatial in- 
terpolation and linear temporal interpolation to get the 
field at the particle position. Initially we release 1.6 mil- 
lion test electrons which have a Maxwellian distribution 
with T e = 2.0 x 10 6 K in the upstream frame. The test 
electrons are released upstream of the shock at the time 
f2 c jt = 30 when the shock is fully formed and far from the 
boundaries. They are released uniformly in space over all 
values of z and between x — 350c/ uj p i and x = 450c/w P i. 
We implement the so-called Bulirsh-S tor method to inte - 
grate the trajectories of the electrons (jPress "eTaI1[T986l) . 
The method is highly accurate and fast when the fields 



are smooth on the scale of particle gyroradius. The al- 
gorithm uses an adaptive method based on the evalua- 
tion of the local truncation error, which is essential to 
reduce numerical error when electrons experience rapid 
field variations. The time step is allowed to vary be- 
tween 10 -3 and 0.257" 1 and the ratio il ce /tt C i is taken 
to be the realistic value 1836. The boundary condition 
in the z-direction is taken to be periodic. The simula- 
tion domain in the x-direction is large enough so that 
no test-particle electrons escape from the system. The 
readers are referred to earlier works using this method for 
additional nu merical details (jKrauss-Varban et alJll989l 
lBurgessll2006ft . Strictly speaking, this test-particle simu- 
lation is only valid when the influence of the accelerated 
electron to the background fluid is negligible. However, 
in the end of the simulation, the electron distribution 
considerably departs from the initial Maxwellian due to 
the energization process at the shock front. The effect of 
the non-thermal electrons on the long-term evolution of 
the termination shock, which is not included in our cal- 
culation, may become significant. We also note that in 
our two-dimensional simulation (or any simulation with 
at least one ignorable coordinate) , the moti on of charged 
parti cle is tied to their original field lines (jJokipii et al.l 
1 1993ft . The effect of particles moving off their field lines 
of force could be important to this process, and requires 
fully three-dimensional simulations, which is beyond our 
available computing resources. 

3. SIMULATION RESULTS 

The shocks in our study have low Mach numbers and 
high upstream inflow speed (corresponding to high speed 
jet from the reconnection outflow region). The plasma 
in the solar corona is hot and has a low beta (strong 
magnetic field) , which is considerably different from that 
in interplanetary space. We focus on the modification 
of the shock surface by upstream Alfvenic fluctuations 
and its effect on the acceleration of electrons. Table Q] 
lists some key parameters for all the simulation runs in- 
cluding the size of the simulation domain, the variance 
of the injected magnetic fluctuation, and the fraction of 
electrons whose energy is more than 15keV at the end 
of simulation (also see Figure 5). For runs 1-4 we con- 
sider the effect of different variances of magnetic turbu- 
lence. The turbulence variances range from 0.0 to 0.3 
and the size of the simulation domain is L x X L z = 
500c/ojpi x 400 c/uj P i (1. 27 km x 1.02km) for these four 
cases. For runs 5-7, the magnetic variances are the same 
as runs 2-4, but the size of the simulation box is changed 
to L x x L z = 500c/ujpi x 800c/wjri(1.27fcTn x 2.03£;m) to 
examine the effect of changing the coherence length of 
the turbulence which, in our case, is governed by the 
size of the box. In flare environment, strong large-scale 
Alfvenic magnetic fluctuations can be triggered by re- 
connection, and cascade to small scales. This process is 
usually assumed to be the source of magneti c turbulence 
required in many accelerati on models fe.g.. IMiller et all 
[l99llPetrosian fc Liul[200l . Since the size of our largest 
simulation domain is still much smaller than the observed 
hard X-ray emission region (L ~ 10 3 /cto), we do not 
consider a realistic geometry of flare termination shock 
but approximate it locally as a perpendicular shock that 
propagates into a plasma containing magnetic fluctua- 
tions. 
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Figure [2] shows (a) the z-component of the magnetic 
field and (b) the ion number density n p , in a color-coded 
representation from run 3 at f2 C it = 110.0. The magnetic 
field and plasma density have been normalized using the 
average upstream magnetic field B w and in-flow den- 
sity no- The averaged Alfven Mach number in the shock 
frame is about 2.0 and the average compressi on ratio is 
about 2 .1. As noted in the earlier works of iGiacalonel 
(|2005bD and lGuo fc Giacalone! (|2010l) . the shock surface 
becomes distorted due to the interaction between the 
shock front and the upstream turbulence. Meandering 
magnetic field lines cross the shock front at various lo- 
cations along the shock, which allows the electrons to 
cross and/or get reflected at the shock front multiple 
times. The shock-front rippling has also been shown to 
contribute to par ticle acceleration by mir roring electrons 
between ripples (|Guo fe Giacalonell2010D . 

Figure [3] presents the energy spectra, dJ/dE vs. E, 
of electrons. The energy spectra are normalized using 
N e vthe{u pi I c) 2 / keV , where N e is the total number of 
electrons used in the simulations and Vthe is the initial 
electron thermal speed. The green solid line shows the 
initial distribution of thermal electrons in the upstream 
region. The black solid line displays the energy distribu- 
tion for all the electrons downstream of the shock at the 
end of simulation (f2 C ji = 130.0) for run 1. In the case 
of no pre-existing fluctuation the electron energization 
is primarily due to heating and nonthermal acceleration 
within the shock layer. The resulting electron spectrum 
is broadened compared to the distribution incident on 
the shock and has a power-law-like tail extending to high 
energy. This is similar to previous work and confirms 
the effect of ion-scale ripp les on the elect ron energization 
at the shock transition ((Burgess 200lf). The effective 
electron kinetic temperature jump across the shock in- 
cluding the supcrthcrmal distribution is about 6 times 
that of the upstream temperature. The electron temper- 
ature jump is about 40% of the proton temperature jump 
across the shock layer in our hybrid simulation. This is 
consistent with the theoretical prediction that the heat- 
ing of electrons in fast shocks is less than that of ions 
(|Goodrich fc Scudderi [l98i IScudderl fl99l and the ob- 
servational constraints from measureme nts at planetary 
bow shocks and interplan etary shocks (jThomsen et al.l 
[19871 ISchwartz et al.lfl98ll . 

In the following we focus on the nonthermal acceler- 
ation of electrons at shocks after considering the pre- 
existing magnetic fluctuations. The blue solid, dot and 
dashed lines in Figure [3] show energy distribution for all 
the electrons in the downstream region at the end of the 
simulation (fl C it = 130.0) for runs 2-4, respectively. At 
this time most of electrons are downstream of the shock 
so the energy spectra do not evolve any longer. It can 
be seen that the electrons are accelerated to high en- 
ergy when the upstream magnetic turbulence is included. 
For higher variance of magnetic turbulence, there are 
more particles accelerated to high energy. For run 4 
(SB 2 /Bq = 0.3), 9.8% of electrons are accelerated to 
more than 15keV by the end of the simulation. The effi- 
cient electron acceleration can be understood as stronger 
magnetic turbulence allows increased field-line meander- 
ing, and the electrons move along field lines of force that 
intersect the shock at several places which therefore al- 
lows for efficient energy gain at the shock front. Com- 



pared with earlier work (|Guo fc Giac alone 201(3), we find 
the maximum energy is proportional to the square of up- 
stream inflow speed in shock frame. This is consistent 
with diffusive shock acceleration in which the energy gain 
is proportional to V 2 h , where V s h is shock speed in up- 
stream frame. 

In Figure 2] we examine the effect of changing the co- 
herence length of the magnetic turbulence and focus on 
the high energy part of the energy spectra. It shows re- 
sults from runs 5-7 (red lines, L z — 800c/ uj p i) along with 
corresponding runs 2-4 (blue lines, L z = 400c/w P i). It 
is shown that for larger coherence length, the electrons 
reach higher energy and the spectral slope is flatter. The 
more efficient acceleration in runs 5-7 can be understood 
as the larger simulation domain in the direction of mag- 
netic field allows more field line wandering normal to 
the shock (AA 2 cx AZ, where AA is the field-line ran- 
dom walk normal to the averaged magnetic field and AZ 
is distance along the field) therefore the electrons move 
across the shock more easily. This dependence shows 
that long-wavelength fluctuations are important to ac- 
celerate electrons to high energy. 

Figure 5 shows the relation between the amplitude of 
magnetic turbulence SB 2 / B 2 injected in hybrid simula- 
tion and the percentage of electrons eventually acceler- 
ated to more than 15 keV. The triangles represent the 
case with L z = 400c/w P i whereas the squares represent 
the case with L z = 800c/ uj p i. This result shows that for 
larger amplitude of magnetic turbulence, more electrons 
get accelerated to high energies. Once SB 2 /B 2 reaches 
0.3 or higher, about 10% of electrons are eventually ac- 
celerated to more than 15keV. For the cases in which 
L z = 800c/o;pi, the simulations generally give slightly 
more electrons accelerated to high energies. 

We also analyze the acceleration of protons, which 
are treated self-consistcntly in this problem (i.e., they 
arc included in the hybrid simulation). Figure 6 shows 
the downstream energy spectra of protons in shock 
frame at fl C it = 130.0 for runs 2-7, normalized using 
NpVthp(w P i I c) 2 1 keV , where N p is the total number of 
protons used to plot the spectra and Vthp is the initial 
proton thermal speed. Similar to Figure El the results 
from runs 5-7 are represented by red lines and the re- 
sults from runs 2-4 are displayed using blue lines. The 
accelerated protons show a similar dependence on turbu- 
lence variance and coherence length to that of electrons. 
This dependence has been found previously in the case of 
higher Mach num ber shock and larger coherence length 
(|Giacalond l2005bT ). These results show both electrons 
and protons are efficiently accelerated. However, for the 
parameters we use, the slopes of the energy spectra of 
protons are considerably steeper than that of the spec- 
tra of electrons. This is probably due to the fact that 
fast moving electrons can interact with shock front more 
times than protons in the given time, which allows a 
more efficient acceleration. Protons may need more time 
and l arger spatial scal es to reach an efficien t accelera- 
tion (IGiacalonel I2005af) . Also, as shown in (j Giacalone! 
l2005bF 7 the effect of pre-existing fluctuations tend to be 
prominent for ions in the case of large coherence length 
(~ 4000c/w P i). Exploring the relative acceleration be- 
tween electrons and ions for different parameters will be 
a subject of future work. 
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4. DISCUSSION AND CONCLUSIONS 

Understanding particle acceleration in solar flares is 
a challenge since only remote observations are available 
and it is hard to identify the main mechanism. While it 
is commonly thought that magnetic reconnection drives 
the energy release, the detailed physical process involved 
in accelerating the electrons and ions is still not clear. In 
this paper we studied electron acceleration at flare ter- 
mination shocks whose existence in the flaring region has 
been predicted by numerical simulations and flare mod- 
els. We find that electrons are rapidly and efficiently 
accelerated at such shocks in the presence of pre-existing 
magnetic fluctuations wit h parameters s i milar to that in- 
ferred by previous work (jForbesi If 9861 : iWorkman et al.l 
1 2 1 If ) . The electrons are accelerated to a few MeV in 
100 ion cyclotron periods (of the order of a millisecond) 
and more than 10% of thermal electrons are accelerated 
to over 15/ceU given a sufficiently strong magnetic tur- 
bulence (SB 2 /Bq ~ 0.3). We also show electron accel- 
eration is more efficient for larger turbulence variance 
SB 2 / Bq and/or a larger turbulence coherence length L c . 
Both of these indicate that large-scale field-line mean- 
dering plays an essential role in accelerating electrons at 
shock fronts. Our simulations suggest that when mag- 
netic turbulence is present the flare termination shock 
could accelerate electrons to much h igher energies tha n 
simple shock-drift acceleration (e.g.. iMann et al1l2009f ). 
We note that a similar mechanism has been shown to 



efficiently accelerate ions and has similar dependence on 
the turbulence properties. This correlation between ions 
and electrons is actuall y commonly observed in solar en- 
ergetic particle events (Oliver 2000). For the parameters 
we use in our simulations, the accelerated protons have 
energy spectra steeper than that of electrons. This is dif- 
ferent from the previ ous results for parameters similar to 
interp lanetary space (|Giacalonel l2005b; G uo fc Giacalonel 
120101) . We note that these results are carried out for en- 
ergies lower than the injection energy for diffusive shock 
acceleration and can be variable and depend on param- 
eters such as the property of magnetic turbulence. 

We also note for the situation we study, the resulting 
distribution of electrons is non-Maxwellian. The struc- 
ture of collisionlcss shock may be considerably modified 
by the accelerated particles. While this effect is not con- 
sidered in our test particle simulations, it may be impor- 
tant. A full particle simulation with similar parameters 
will also be useful to study the energy partition between 
electron and protons during the acceleration at the flare 
termination shock. While our results are qualitatively 
consistent with the observed hard X-ray emissions in so- 
lar flares, other plasma effects in the flare region may 
need to be considered to directly compare with the ob- 
servations in the future. 
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Run 


T (r 1 ti) -\ y Ti~( r hi) -\ 


SB 2 IB? 


r%(E > isfrpVi 


1 


500 x 400 


0.0 


1.3 


2 


500 x 400 


0.03 


4.5 


3 


500 x 400 


0.1 


8.0 


4 


500 x 400 


0.3 


9.8 


5 


500 x 800 


0.03 


4.9 


6 


500 x 800 


0.1 


8.9 


7 


500 x 800 


0.3 


11.9 



TABLE 1 

Some parameters for different simulation runs. The size of the simulation domain, the variance of injected magnetic 
fluctuation, and the fraction of electrons whose energy is more than ibkev at the end of simulation. 




Fig. 1. — A cartoon illustration of geometry of flare termination shock. 
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Fig. 2. — The color-coded images of (a) magnetic field B z , (b) ion density n p for run 3 at fl c it = 110.0 
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Fig. 3. — The energy spectra of electrons at the end of the simulation (fl^t = 130.0). The energy spectra are normalized using 
NeVthei^pi I c ) 2 1 keV , where N e is the total number of electrons used in the simulations and v^e is the initial electron thermal speed. 
The green solid line shows the initial distribution of thermal electrons in the upstream region. The black solid line displays the energy 
distribution for all the electrons in downstream region at the end of simulation for run 1. The blue solid, dot and dashed lines represent 
results from runs 2, 3, and 4, respectively. 
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Fig. 4. — The energy spectra of electrons at the end of the simulation (flit = 130.0). The energy spectra are normalized using 
N e v t f le (id p i/c) 2 /keV , where N e is the total number of electrons used in the simulations and v the is the initial electron thermal speed. 
The red solid, dot and dashed lines represent the energy distributions for all the electrons in downstream region at the end of simulation 
for runs 5, 6, and 7 respectively. The blue solid, dot and dashed lines represent results from runs 2, 3, and 4, respectively. 
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Fig. 5. — The relation between the turbulence amplitude SB 2 / B' 2 , injected in hybrid simulation and the percentage of electrons eventually 
accelerated to more than 15 keV 
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Fig. 6. — The energy spectra of protons downstream of the shock at the end of the simulation {Q.it = 130.0), normalized using 
NpVth-pi^pi / c) 2 / keV , where N p is the total number of protons used to plot the spectra and vthp is the initial proton thermal speed. 
The red solid, dot and dashed lines represent the energy spectra for protons in downstream region at the end of simulation for runs 5, 6, 
and 7 respectively. The blue solid, dot and dashed lines represent results from runs 2, 3, and 4, respectively. 



